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We present a scheme for rapidly entangling matter qubits in order to create graph states for one-way quantum 
computing. The qubits can be simple 3-level systems in separate cavities. Coupling involves only local fields 
and a static (unswitched) linear optics network. Fusion of graph state sections occurs with, in principle, zero 
probability of damaging the nascent graph state. We avoid the finite thresholds of other schemes by operating 
on two entangled pairs, so that each generates exactly one photon. We do not require the relatively slow single 
qubit local flips to be applied during the growth phase: growth of the graph state can then become a purely 
optical process. The scheme naturally generates graph states with vertices of high degree and so is easily able to 
construct minimal graph states, with consequent resource savings. The most efficient approach will be to create 
new graph state edges even as qubits elsewhere are measured, in a 'just in time' approach. An error analysis 
indicates that the scheme is relatively robust against imperfections in the apparatus. 
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I. INTRODUCTION 

Despite significant and exciting experimental progress in 
recent years, the physical realization of a full-scale quan- 
tum computer (QC) remains a tremedous challenge [1]. In 
many systems excellent single qubits have already been real- 
ized (notably, of ions in a trap [2, 3], NV centres in diamond 
[4, 5], etc). However, few systems have demonstrated con- 
trolled qubit-qubit coupling between pairs taken from more 
than four qubits, and achieving the necessary exquisite con- 
trol remains highly problematic. In general it is difficult to 
simultaneously satisfy the two key requirements of coupling 
diffent subsystems in a controlled manner, while at the same 
time shielding the system from its environment [1]. In the 
majority of QC schemes, some direct physical interaction is 
supposed to generate the two-qubit operations (e.g., phonon 
modes among trapped ions [2, 3], Forster interactions between 
excitions in semiconduction quantum dots [6], etc). Thus one 
calls for the qubits to strongly interact with selective parts of 
their environment (namely, other qubits and the control mech- 
anisms) while avoiding interactions the rest of the environ- 
ment to a near perfect degree. This is obviously an challeng- 
ing prescription. 

As an alternative to employing a direct physical interaction 
between qubits, one can exploit the entangling power of mea- 
surements. A suitable measurement, at least for certain out- 
comes, will have the effect of projecting previously separate 
qubits into a highly entangled state. This idea has been ex- 
plored as a route to QC using photon qubits in a linear optical 
apparatus. Measurement-based gates have indeed been shown 
to be sufficient for universal gate-based quantum computation 
[7]. However, in order to achieve each logical gate with high 
probability, one must prepare and then consume large auxil- 
iary resources. This necessity is essentially due to the small 
probability of success of the elementary quantum gates based 
on auxiliary systems and measurements [8]. 

One way to reduce this overhead is to exploit the idea of 
one-way computing [9, 10]. In this approach one would pre- 



pare a certain multi-qubit entangled state, a cluster [9] or a 
graph state [10, 11], prior to the computation. This state has 
the property that the computation can then proceed purely by 
single-qubit measurement - essentially consuming the graph 
entanglement as a resource. Recently there has been a suc- 
cessful proof-of -principle experiment realising a 4-qubit clus- 
ter state [12]. A key advantage of the one-way computing 
strategy is that it introduces a degree of separation between the 
act of creating entanglement and the act of executing the com- 
putation. Thus we need not expend the effort needed to ensure 
that each entangling operation succeeds with high probability 
- we can tolerate failures during the growth process simply 
by rebuilding the affected graph section, provided of course 
that failures are heralded. Indeed, in this spirit various re- 
cent schemes [13, 14, 15] have shown how to take gate oper- 
ations that are fundamentally non-deterministic, and use them 
to construct an such an entangled resource state with certainty. 

One particularly attractive possibility is to use matter 
qubits, with the obvious benefits that they are static and po- 
tentially long lived, together with an optical coupling mech- 
anism that creates suitable entanglement. Based on earlier 
schemes that allow for generating entanglement or realizing 
quantum gates in matter qubits using flying optical qubits 
[16, 17, 18, 19, 20, 21, 22, 23], two recent publications 
[24, 25] in particular have explored precisely this possibil- 
ity. The matter qubits can be completely separate, for exam- 
ple each within its own cavity apparatus, providing that suit- 
able optical channels connect them to a mutual measurement 
apparatus. The simplest scheme is that of Barrett and Kok 
(BK) [24], where one requires only a single beam splitter and 
two detectors in order to couple pairs of qubits. The elegant 
BK approach however suffers from the constraint that, even 
with ideal apparatus, the entangling operation must fail with 
a probability of p = 1/2. Failures damage the nascent graph 
state, but because the failure is flagged, or 'heralded' , the dam- 
aged parts can be removed and the growth can continue. Nev- 
ertheless, the high rate of destructive failures introduces a con- 
siderable overhead [26], especially with certain types of target 
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graph topology. The scheme due to Lim, Beige, and Kwek 
(LBK) [25] introduces the idea of 'repeat until success' entan- 
glement, meaning that while failures still occur with probabil- 
ity p > 1/2, these failures are essentially passive and one can 
simply try again. Thus one can construct graph states with a 
lower overhead, in terms of number of entangling operations, 
and any topology can be directly implemented. However, the 
cost for this advance is that the underlying coupling process 
is more complex: each matter qubit gives rise to a superpo- 
sition of an 'early' and a iate' photon in time-bin encoding, 
which must subsequently enter a beam splitter apparatus si- 
multaneously. This appears to be more challenging relative to 
the simpler BK scheme, so that it is an open question which 
scheme is the more practical. 

Here our goal is to unite the more desirable features of both 
these schemes, in particular the simple static optical appara- 
tus of the BK scheme and the non-destructive 'repeat-until- 
success' aspect of the LBK approach. Moreover we introduce 
a vital feature which neither of these approaches possess: we 
demonstrate a graph growth mechanism which does not re- 
quire local unitary operations (e.g., flips) to be performed on 
the matter qubits during the growth process. The growth then 
becomes purely a sequence of optical excitations, with a cor- 
responding significant increase in speed and considerable re- 
duction in complexity. 

We intend that the present paper will form a self contained 
overview of the entire paradigm that we are advocating, and to 
this end we include compact analysis of the relevant proper- 
ties of graph states. We make use of the idea of a minimal 
graph state (MGS), and make a comparison with the more 
limited 'cluster states' which results when the geometry of 
physical qubits and their neighbors are fixed by experimental 
constraints. We conclude that there are dramatic savings, in 
terms of qubits and entangling operations, when one adopts 
an architecture that can build an MGS directly. 

II. GRAPH STATES AND CLUSTER STATES 

Graph states [10, 11, 27] are multi-qubit entangled states, 
which can be conceived as having been entangled according to 
certain pattern of two-qubit phase gates. Formally, this pattern 
is specified by the adjacency matrix of an (undirected simple) 
graph G(V, E), where V denotes a set of n vertices associated 
with the qubits, and edge set E reflecting the phase gates (see 
Fig. 1 (d) for example). The graph state of the empty graph 
has the state vector |*) = \+}® n = ((|0) + \1})/V2)® n . 
The state vector of a graph state including edges can then be 
written as 

\g) = n p( a < b )|+)®", (i) 

{a,b)eE 

with p( a ' h ) corresponding to a phase gate p( a > 6 ) = (l+cri°' ) + 

)/2 between qubits labeled a and b, expressed 
in terms of Pauli operators. Such graph states are stabilizer 
states [28], and in turn, every stabilizer state of n qubits is 
locally equivalent to a graph state [29, 30]. 



A cluster state (CS) [9] is a particular graph state: it is one 
with an underlying cubic lattice of one, two or three dimen- 
sions (see Fig. 2 (a) for example). A cluster state of more than 
one dimension has the remarkable property that it forms a uni- 
versal resource for measurement-based one-way computing: 
having created this state, the actual computation is executed 
simply by making local measurements [9, 10]. It is universal 
in the sense that the procedure amounts to effectively imple- 
menting an arbitrary unitary on the input qubits. 

However, the measurements performed in order to imple- 
ment some chosen algorithm will include two classes which it 
is important to distinguish [9, 10]. The first are the Pauli mea- 
surements, which we can denote as measurements along the 
X, Y, or Z axis. Each such measurement maps a graph state 
onto another graph state for all outcomes. For example the 
Z measurement effectively deletes the measured qubit (node) 
and its associated edges, while X and Y measurements alter 
the graph according to the rules given in Ref . [11]. These mea- 
surements correspond to the Clifford-part of the computation, 
and the resulting map on the level of states can always effi- 
ciently be determined on a classical computer [10]. Having 
performed all the prescribed Pauli measurements on a cluster, 
we are left with a minimal graph state (MGS) which is the 
graph containing the smallest number qubits that is capable 
of realizing our desired algorithm. The remaining measure- 
ments are of the second class: von Neumann measurements 
in tilted bases. Such measurements take the system out of 
the graph state prescription and generally cannot be efficiently 
simulated on a classical computer. In a sense one can think of 
the Pauli measurements as simply customizing the (initially 
universal) cluster state into the form that will implement our 
chosen algorithm, while the more general tilted measurements 
actually execute the algorithm. 

Many physical systems that can generate graph states are 
in fact limited to cluster state generation, because the physical 
qubit interactions are limited to some kind of nearest-neighbor 
(or at any rate, local) form. This applies to implementations 
in electron spin lattices and optical lattices. However, we are 
under no such constraint since the physical qubits have no de- 
fined geometry [31]. Instead, we can directly 'grow' an ar- 
bitrary graph, and hence we may prepare the graph state that 
forms the specific resource for a given quantum algorithm. 
We would therefore seek to directly build a MGS, shortcutting 
the creation of the cluster state with its redundant universality. 
This proves to have dramatic advantages in terms of the num- 
ber of entanglement operations and qubits needed. In general 
one finds that a MGS will often exhibit a high vertex degree, 
and will be contain significantly fewer qubits compared to the 
graph state that is obtained from a cluster state after measure- 
ments along the Z basis, essentially merely removing qubits 
(typically up to an order of magitude). Explicit examples are 
described later. 



III. DESCRIPTION OF THE PHYSICAL SCHEME 

In this section we describe the physical requirements and 
processes involved in implementing our proposal. We start by 



describing the elementary physical systems required. We then 
outline the action of the beam splitter network, both on simple 
product states and more importantly when fusing graph states 
together. Finally we then show how to do this without any 
single qubit unitaries, and make some concluding comments. 



A. Physical Components 

In Fig. 1 (a) we indicate the basic energy level scheme that 
our matter qubits should incorporate. Obviously, a real quan- 
tum system may have additional levels, but providing these 
3-level dynamics are incorporated, then such a system is suit- 
able. Candidates include an actual atom or ion in a trap, of 
course, but we may also consider any optically active solid 
state structure with a discrete spectrum, such as a quantum 
dot, or NV-diamond centre [40]. The ground states, labeled 
with state vectors |0), |1), are the qubit basis states. The 
third level, labeled |e), provides a mechanism for produc- 
ing a photon from the atom, conditional on the atom being 
in state labeled with state vector That is, there is an ex- 
ternally driven transition from |1) — > |e) by using a 7r-pulse, 
followed by the optical relaxation |e) — > |1) which emits a 
single photon into the cavity mode, and eventually 'leaks' out 
to an external optical network. In a single-mode description, 
this is characterized by a coupling strength O of the Jaynes- 
Cummings coupling between the transition |1) — * |e) and the 
cavity mode, with decay rate T of the cavity mode. The sys- 
tem is continuously observed via the photon detectors placed 
behind the four beam splitter array (4BS). We note that, as re- 
marked in Ref. [25], if we have a fourth level accessible from 
|e), then we can potentially create a photon directly in the 
cavity without significantly populating |e) - this may be ad- 
vantageous in avoiding dephasing. The state labeled |e) is also 
exploited when we wish to make a measurement - continuous 
illumination by a laser adjusted to the transition energy will 
result in a fluorescence conditional on the qubit state. This is 
a Z measurement; measurement in the other directions is ac- 
complished by an appropriate local rotation followed by this 
fluorescence. 

The elementary multi-qubit operation in our proposal is 
based on a four-port beam splitter (4BS), which is a composite 
of four ordinary beam splitters, arranged so that every input 
crosses every other, and finally incident on photon counters. 
Two types of 4BS will be employed: one without additional 
phase shifts, the basic network, and one which includes a cer- 
tain phase-shift corresponding to a factor of e i7r / 2 , the shifter 
network, for producing certain important cluster states. The 
action of the 4BS is essentially to 'erase' information about 
which cavity a photon originated from, so that a given de- 
tector cannot differentiate between matter qubits. Ideally the 
frequencies of the modes are identical, and other sources of 
mode matching problems are to a high extent eliminated - we 
analyze the effects of realistic imperfections in Section VI. 
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FIG. 1 : (a) The energy level scheme for a matter qubit. (b) Our ap- 
paratus: matter qubits j — 1, 4 (lower) emit photons to detectors 
k = 1, 4 (upper) via beamsplitters. We consider two variants of 
the device: one with a phase shifter as marked, one without, (c) The 
effect of excitation and measurement when four product state qubits 
are employed - there are various possibly results depending on the 
number and pattern of detected photons, (d) A graph state in which a 
'leaf qubit, marked b, is attached by only one edge, (e) The effect of 
applying our protocol (with shifter) to two EPR pairs, (f) The effect 
of applying our protocol (without shifter) to two arbitrary graphs as 
shown. 



B. Action on Product States 

The analysis to determine the specific projections is 
straightforward. The most simple interesting case we examine 
is that of inputting four qubits in the product state correspond- 
ing to |+)® 4 . The action of the optical excitation |1) — > |e) 
applied to all qubits, followed by the emission into their lo- 
cal cavities, then results in an equal superposition of all basis 
vectors containing all binary words, 

\<P) = \ E \i,j,k,lKc\y(4) j (4) k (cl) l \0). (2) 

where the annihilation operators of the respective cavity 
modes are denoted by c%, C4. 

As the photons propagate through a beam splitter into new 
modes we employ mappings such as (dj, d\) T = B(cj, c\.) T ', 
j,k = 1, 4, where phases of the transmitted and reflected 



4 



mode are chosen such that B is given by 

Thus the network corresponds to a unitary manipulation of the 
photon states. In this way we eventually obtain the final gen- 
eration of operators representing photons in the modes upon 
which the detectors act. Thus we determine the state that re- 
sults from a particular detector reading. In fact a number of 
states can occur, including states that are locally equivalent to 
linear 4-qubit graphs and 3-nodes, as shown in Fig. 1 (c). 



C. Fusion of Graph States 

The ability to generate such graphs directly from product 
states appears to be a promising characteristic. However, 
to properly differentiate the possible outcomes one would 
require either resolving photon detectors (capable of distin- 
guishing a single photon from a pair, etc) or else one would 
need to resort to a lengthy asymptotic variant of the 'double 
heralding' idea in Ref. [24]. These undesirable features result 
from the fact that the systems state, prior to photon measure- 
ment, was not an eigenstate of the total photon number oper- 
ator: there are elements in the superposition corresponding to 
0, 1, 4 photons. To avoid the problem we must contrive to 
introduce a known number of photons. This is the same issue 
faced in Ref. [25], where the authors suggest resorting time- 
bin approach and local flips in order to guarantee that each 
matter qubit ultimately generates one photon. 

We take a different route, based on the idea of fusing graphs 
together. We will find that we can regard EPR pairs as a kind 
of raw ingredient from which graphs of arbitrary complex- 
ity can be grown deterministically. Recall that an EPR pair 
with state vector |EPR) = (|0, 1) - |1, 0))/\/2 is already LU- 
equivalent to the simplest non-trivial graph state, the one con- 
sisting of two vertices connected by an edge. We use the term 
LU- equivalent to mean, equivalent up to local unitary oper- 
ations on individual qubits. Our fusion process exploits ex- 
isting entanglement within the graphs sections: certain vertex 
pairs within a graph state can be locally rotated to the sub- 
space span{|0, 1),|1,0)} - two such pairs then generate pre- 
cisely two photons. 

Suppose that a 'leaf node exists, i.e., a certain vertex (as- 
sociated with qubit a) is attached to only one other vertex (as- 
sociated with qubit b) of the graph. This is shown in Fig. 1 
(d). Then the state vector \^g) of the entire graph state is of 
the following form 

|* G ) = (10, 0) + [1, 0) + |0, l)P - |1, 1)P) |V) (4) 

where the vectors \a, b) correspond to the qubits a, b, and \ip) 
refers to the external, arbitrarily connected part of the graph 
state, shown inside green bubbles in Fig. 1 (d). We define 

p= n ^< ^ 

i£N b \{a} 



so with index i running over the neighbors Nb of qubit b lying 
within \ip). This state vector is in turn equivalent, up to a 

unitary rotation ( 1 — ia^ ) / \/2 on a, to 

|* s ) = (|l,0) + |0,l)i 5 )|V) (6) 

Having made this transformation, we know that our qubit pair 
labeled a and b will emit precisely one photon. If we similarly 
prepare a second pair of qubits, associated with a different 
graph (or, a different part of the same graph) then the two pairs 
of qubits can be employed in our4BS device and will generate 
precisely two photons. As indicated in Fig. 1, we have con- 
sidered two variants of the beam splitter network: one with 
and one without a phase shifter. A simple analysis determines 
that in both cases there are two possible classes of outcome. 
The two photons may arrive in a single detector, in which case 
the effect is simply equivalent to applying local phase gates to 
the matter qubits. Alternatively, the two photons may enter 
different detectors, in which case the two pieces of graph are 
fused together in a fashion we specify presently. Identifying 
the various outcomes does not require counting photon detec- 
tors, since we know there are two photons in total (but such 
detectors may be useful in fighting errors when taking imper- 
fections into account, see Section VI). 

The two classes of outcome are equally probable. In the 
case of the former one can try again without pausing to cor- 
rect the local phases, which can be fixed after the eventual 
successful fusion. The average number of attempts required 
is two. This is then a 'repeat-until-success' scenario equiva- 
lent to the one first observed in Ref. [25]. The particular form 
for our fused graph depends on whether the phase shifter was 
employed. If we do employ the shifter, and supposing that we 
input two EPR pairs, then the resulting state is LU-equivalent 
to a linear four-qubit graph as shown in Fig. 1 (e). If instead 
we use our 4BS without the phase shifter, i.e. the basic net- 
work, we can couple arbitrary graph fragments according to 
the rule shown in Fig. 1 (f). For the example of joining two 
EPR states, we show the outcome states conditioned on which 
detectors click, and their probabilities in Table I. These obser- 
vations lead us to regard EPR pairs as the basic resource for 
graph growth: EPR pairs can be generated easily by a single 
beam splitter using the BK scheme, or equivalently we have 
observed that they can be obtained from our 4BS network by 
choosing to excite just two of the four qubits. The combina- 
tion of the two coupling processes Fig. 1 (e) and Fig. 1 (f) then 
allows graphs of arbitrary complexity to be built. Recall that 
a graph can be 'pruned', i.e., qubit nodes can simply be re- 
moved, by making a Z axis measurement, while other useful 
transforms result from X or Y measurements [11]. Indeed, 
a recent preprint [41] makes ingenious use of measurements 
on leaf structures, reminiscent of those occurring at the fusion 
point in Fig. 1(f), for qubit loss tolerance in graph states. 



D. Growing a Graph Without Employing Local Gates 

We have seen that we can grow graphs by transforming se- 
lected qubits to an EPR-type basis prior to fusion, and then 
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1>(fci;fc 2 ));Prob | 12 3 4 

~T~ ~ |a);I 0;0 |c>; ^ \d); i 

2 0;0 |a);| |d>; £ |c>; £ 

3 |c);£ |d);& |6);i 0; 

4 |c);i 0;0 |6>;± 

TABLE I: States vectors resulting from clicks in detectors fci and fc2, 
and their probabilites, for the basic beam splitter network acting on 
the state |EPRi, 2 )|EPR 3 ,4), where \a) = (|0. 1) +i\l, 0})(|0, 1) - 
i|l,0»/2, 1 6) = (|0,1) - <|1,0»(|0,1) + <|l,0»/2, |c) = 
(|0, 1, 0, 1) - jl, 0, 1, 0))/y/2, \d) = (|0, 1, 1, 0) + |1, 0, 0, 1})/V2. 
State vectors |c) and \d) are LU-equivalent to a graph state vector in 
which a central vertex radiates three 'leaf vertices. 

applying additional LU operations to transform the 'raw' re- 
sultant state back to a graph state. But, can we avoid these lo- 
cal transformations? It is evidently necessary to employ single 
qubit rotations at two stages: the very beginning in the entan- 
gling procedure, where we must take 'fresh' qubits and pre- 
pare them in |+) in order to synthesize the EPR pairs which 
we regard as our basic ingredient, and the very end where we 
will wish to rotate qubits prior to our fluorescence measure- 
ment, in order to synthesize measurements along some gen- 
eral axis. Remarkably, we can in fact omit the numerous local 
rotations during graph state growth. We find that, within a 
light constraint on the growth process, we can ensure the state 
remains LU-equivalent to a graph state at each growth step. 

To see that this is possible consider the following argument. 
Suppose that we have some multi-qubit state vector \ ty) which 
meets the following two conditions: 

(i) The state vector |^) is equivalent up to local unitary 
operations to \^g) corresponding to a graph G of the form 
as in Fig. 1 (d). 

(ii) Regarding the pair of qubits labeled a and b, 

|(0,1|*)| = |(1,0|*)|, K0,0|*>| = |<1,1|*>|=0. (7) 

From (i) and recalling Eqn. (6) we know that \ty) is LU- 
equivalent to a state vector | 'J e) = 0) + |0, 1}P)\i/j) since 
that is itself LU-equivalent to the graph state vector corre- 
sponding to G. The additional constraint (ii) implies that our 
state vector can be written as 

|tf) = (|l,0)£/+|0,l)£/)hi>) (8) 

where U is a product of local unitaries acting on i.e., act- 
ing in the Hilbert spaces of the qubits other than those labeled 
a and b, and U = exp(icf))UP with (f> 6 [0, 2tt) an arbi- 
trary phase. Now let us apply the 4BS process of Fig. 1 (f) 
to a and 6, along with an equivalent pair from some analogous 
graph state (or, another part of the same graph state). On fail- 
ure (with probability 1 /2) the process simply introduces some 
known local phases, which do not alter our prescription. On 
eventual success we generate a state vector 

|*tot) = (| Wot + |X)[7 tot )|V>tot>. (9) 

Here X is a binary word with two zeros, two ones, and X is 
its complement. The vector |\&tot) refers to the entire state 



vector of the fused system, and |V>tot) is the state vector for 
all qubits except the four that coupled via the 4BS. Similarly, 
t/tot is some product of local unitaries on iV'tot), and [/tot = 
exp(i6>)P to t£/tot, where 

i 

with the index running over neighbors of either of the orig- 
inal vertices (but excluding mutual neighbors). The phase 
9 6 [0, 27r) is determined by <\> and its counterpart in the 
second pair, together with phases introduced in any failures 
preceding the successful fusion. This state vector |\l/tot) is in- 
deed LU-equivalent to the desired fused graph state of Fig. 1 
(f). Moreover, because X has two zeros and two ones, if we 
nominate one of those four qubits to be a new 'vertex' qubit 
b, two of the three remaining 'leaf qubits are available to be 
labeled as a to satisfy (i) and (ii). Thus we can go on to per- 
form further fusions using l^tot)- To conclude the argument 
we need only observe that conditions (i) and (ii) are met by 
simple EPR pairs, and by the state resulting from fusing two 
EPR pairs via the process depicted in Fig. 1 (e), i.e. the state 
that is LU-equivalent to a linear four qubit graph state. Thus 
these simple states can act as the initial building blocks as we 
construct a complex graph. 

Then provided we are prepared to measure out one in every 
three of the 'leaf nodes which occur at each fusion point, we 
can grow our entire graph from EPR pairs without the use of 
local unitary operations during the growth process. This con- 
straint is extremely light: we would rarely wish to use all three 
leaves, and in any case the number of leaves can be increased 
by two simply by fusing an EPR pair, which adds one to the 
number of leaves that are eligible in the sense of property (ii). 
Of course, the LU operations needed to map the final state to 
the desired ultimate graph state can subsumed into the rotation 
which in any case precedes measurement. The growth pro- 
cess is therefore entirely one of optical excitation and detec- 
tor monitoring. One can anticipate than in many systems, the 
cost in efficiency arising from following the constraint would 
be vastly outweighed by the increase in growth speed. 



E. Further Remarks on Graph Growth 

There is one additional comment to make regarding the 
speed of our protocol: in the scheme of Ref. [24] it is al- 
ways necessary to wait a period after the initial measurement 
to ensure there are no further photons in the apparatus. The 
fidelity of the entangled states is only high if this wait period 
is long compared to the typical time for a photon to be de- 
tected. This additional waiting time necessary in Ref. [24] is 
not necessarily long, given that photon emission from these 
sources is approximately exponential, governed by the time 
scale 1/Tsiow By contrast, because we contrive to have pre- 
cisely two photons in the apparatus, once we see two detection 
events (either in different detectors, or, given resolving detec- 
tors as discussed later, within one detector) we have no need 
for such a wait. One should hence expect a factor of about 5 in 
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difference concerning the speed of this step in the respective 
schemes. 

In a mature form of the architecture described here, one 
would envisage coupling the many qubits by a form of Al- 
port all-optical router, as used in integrated optics, so that our 
qubits can remain static and we can choose which of them 
will couple by suitably setting the router and optically excit- 
ing only that subset. Devices relevant to this technology have 
already been developed for classical optical communications 
[35]. This would permit direct growth of graphs with an arbi- 
trary topology, and in particular the ability to directly entan- 
gle arbitrary qubits gives a non-local architecture [36, 37, 38]. 
Such non-local architectures may prove to have an advantage 
in quantum fault tolerance and error correction [39]. 



IV. CHARACTERISTICS OF THE DERIVED GRAPH 

STATES 

As discussed above and illustrated in Fig. 1, our protocol 
can generate graphs of arbitrary topology, including nodes 
of high degree. We argued earlier in Section II that, given 
such an architecture, one should aim to create minimal graph 
states (MGS). The advantage in terms of resource consump- 
tion when preparing appropriate MGS compared to standard 
cluster states (CS) can be quite significant. 



A. Minimal Graph States 

In the best known scheme for an n-qubit quantum 
Fourier transformation [10], the number of required qubits is 
CFourier(n) = 8n 2 + 0(n) for a cubic cluster, followed by 
first X, Y, and Z measurements, then tilted measurements. 
The MGS embodies G Fom - isl (n) = (3/2)n 2 + 0(n) qubits. In 
turn, for the quantum adder, 

CW(») =312n + 0(l), G A dd a (n) = 16n + 0(1) (11) 

[10]. Hence, one may gain more than an order of magnitude in 
resource consumption. For the 3-qubit Toffoli gate that we use 
here as our illustration, we have Gxoffoii = 13 versus CVoffoii = 
65, so a factor of 5 difference in the number of qubits. 



B. Generation of Edges 

Concerning the actual preparation of the graph states, we 
emphasise two points: firstly, when introducing edges with 
a physical interaction, one should always prepare the LU- 
equivalent graph state corresponding to the graph with the 
minimal number of edges. Or, more specifically, the graph 
state with the minimal number of edges that is equivalent up 
to local Clifford unitaries [42], which merely amount to a local 
Clifford basis change. This has also been emphasised in Ref. 
[45]. Fortunately, an efficient algorithm is known to check 
full local Clifford equivalence [46]. Any graphs that corre- 
spond to local Clifford equivalent graph states can be related 




FIG. 2: Progressive measurement on qubits (marked red), (a) A fully 
connected cluster state. After Z measurements (b) and X and Y 
measurements we obtain a minimal graph state (c) for the algorithm, 
here a Toffoli gate. In practice one would further reduce the number 
of qubits required by just in time graph creation - any qubit whose set 
of edges are complete can be measured and 'recycled' to the 'front'. 
Side panel (e) indicates the overhead required in the BK scheme in 
making high degree vertices: one measures out a portion of the qubits 
(in the X basis here) in order to alter the graph topology. 



to each other with a successive application of local comple- 
mentations [11, 46]. Also, it is known how many different 
graph states are contained in an equivalence class with respect 
to local Clifford unitaries [29]. 

Secondly, the present scheme seems particularly suitable to 
prepare graph states of graphs involving vertices with a high 
vertex degree in a single step. In a CS after measurement of 
the unused qubits along the Z direction, it suffices to have 
vertices with a maximal vertex degree of 3. This is obviously 
the lowest possible for a graph less trivial than a linear clus- 
ter state. From the example of Fig. 2 one would suspect that 
a typical MGS may need higher degree nodes, and this is an 
important question in considering how they can be efficiently 
constructed. To explore this point we consider the 'maximal 
vertex degree', by which we mean, the highest degree of any 
vertex in the graph. The vertex degree in a MGS can in princi- 
ple take any value. The maximal vertex degree is notably not 
invariant under local Clifford unitaries [42]. To render the no- 
tion of maximal vertex degree meaningful, we have to take its 
minimum value when minimized over all local Clifford uni- 
taries. For a GHZ state of n qubits, it can easily be shown 
that it has a smallest maximal vertex degree of n — 1. So 
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we immediately see that it is meaningful to talk about 'highly 
connected graph states'. For the resource state for the [7,1,3]- 
CSS code as considered in Ref. [11] we find that in the whole 
orbit under local Clifford unitaries the smallest maximal ver- 
tex degree is 6, whereas the largest is 34. In the three-qubit 
quantum Fourier transform the smallest maximal vertex de- 
gree is 4. Thus we see that high degree vertices are indeed 
generally unavoidable in a MGS: any scheme that claims to 
be able to directly and efficiently construct such a state must 
be able to create graph states with vectices of high vertex de- 
gree. The BK scheme, for example, appears somewhat limited 
by the increased difficulty of making the high degree vertices 
associated with graph states; the high rate of destructive fail- 
ures leads one to take an indirect approach as depicted in Fig. 
2(e), with an associated cost in resources. The scheme pre- 
sented here is among the few that generate high degree nodes 
directly [31], as one can quickly see by considering the fusion 
rule depicted in Fig. 1(f). 

C. Cluster States 

Of course, if for some reason we wished to generate a con- 
ventional CS rather than the MGS specific to some given al- 
gorthim, then we can do so efficiently. As an exercise, let us 
conclude this section with a comment on the required number 
of steps in the preparation of a CS with an underlying two- 
dimensional cubic lattice. We will count resources in terms 
of the number of applications of the shifter iVshifter> with two 
EPR pairs fed in in each instance, and of the basic network 
-^Basic- The basic building block can be taken to be a cross 
shape of length 4, requiring four EPR pairs, and the applica- 
tion of two shifter and one basic networks. One row of width 
n can be build using In invokations of the shifter network and 
In uses of the basic network. A two-dimensional cluster state 
onanxn cubic lattice hence requires 

iVshifter(n) - 2n 2 , N B . dsic (n) = 3n 2 - n. (12) 

V. JUST IN TIME GRAPH CREATION 

In a lattice system, one may be well advised to prepare the 
multi-particle cluster state in one step, exploiting a natural 
nearest-neighbor interaction. However, in a scheme as con- 
sidered here, there is no motivation to prepare graph edges far 
in advance of the eventual measurement operations that will 
consume them. One should therefore avoid doing so since 
this gives rise to unnecessary errors due to the graduate degra- 
dation in phase integrity from decoherence. Instead, one can 
introduce new edges and vertices for our MGS shortly before 
it is needed, in a manner analogous to the block-by-block pro- 
cess of Ref. [15] but at a finer scale. By analogy to the term 
used in classical computing, this may be referred to as just 
in time graph state generation [47] (see also Refs. [10, 48]). 
As noted earlier in Section EH, although we may require lo- 
cal unitaries to create the EPR pairs which constitute our 'raw 
ingredient', the remaining steps involved in generating new 
graph structure can take place without such manipulations. 



One can easily confirm that this is possible, even though 
the measurements on earlier parts of the graph are tilted and 
therefore have taken the system to a non-graph state. Con- 
sider a graph state with graph G = (V, E) corresponding to 
the whole computation: let us consider the state vector af- 
ter k measurements on vertices ai,...,afe, forming a vertex 
set T4 C V. The resulting state vector after measurements 
in direction ru, - depending on the measurement outcomes 
Si, ...,Sfc G { — 1,1}, - in this temporal order is given by 
Pk\G) =n- = i(l + (-l) s ^ J -(si,S2,-,^-i)^' ) )/2)|G), 
where a^ a ^ is the vector of Pauli matrices at vertex la- 
beled Oj. Note that the appropriate measurement basis 
fj(si, S2, Sj-i) at step j depends on the earlier measure- 
ment outcomes. Yet, at this point we could have just prepared 

\G k )= n p(°' 6 )|+>®" (O) 

(a,b)e£ fc 

before performing the above measurements, where E). = 
{(a, b) € E : a € Vfc or b £ 14}. Thus the only constraint on 
this just in time approach, is that one should ensure that all 
edges in Ef. are appropriately entangled in step k, see Fig. 2 
(d). 

VI. ERROR ANALYSIS 

A physical implementation of this scheme would be sub- 
ject to a number of possible errors. Our protocol relies on 
the subsystems being identical, so that their outputs are in- 
distinguishable. Thus, mismatching parameters will lead to a 
reduction in performance. Other errors include dephasing of 
the matter qubits, imperfect optical excitation, phase noise (or 
drift) in the optical apparatus, and photon loss, section the ma- 
jority of our analysis will focus on errors due to mismatched 
subsystems; we will comment on the other error sources at the 
end. 

Since the results described here involve the detection of two 
photons arriving from a source, the qualitative effect of errors 
will be similar to the results presented in Ref. [40], and we 
analyse the system using similar methods. For the purpose 
of this analysis, we assume that each atom is a three level 
system, with degenerate ground states, labeled |0) and |1), 
and a level |e) that is optically coupled to |1), with an energy 
hoj e . The cavity is taken to have a frequency uj c = uj e + A, 
where A is nominally zero. The transition |e) *-> |1) couples 
to the cavity mode with a strength SI, and the cavity mode 
decays with a rate T. Thus, we consider here imperfections in 
Aj, ilj and Tj for each atom-cavity subsystem, j = 1, 4. 
A comparable analysis was performed in Ref. [24], so that 
we can compare that two-qubit scheme with the present four- 
qubit protocol. Remarkable, we find that the sensitivity to 
defects in the apparatus is essentially the same. 

A. Continuous Measurement Analysis 

In the following we describe the dynamics of a three-level 
atom in a leaky cavity in the Schrodinger picture, continu- 




FIG. 3: Probability density function, p(t; At), for measuring the first 
photon at time t followed by the second after a delay of At. We have 
used SI = 1, T = 4 and A = for all subsystems. 

ously monitored by a photodetector. Its stochastic dynam- 
ics under continuous measurement can be described using a 
quantum-jump approach, leading to a piecewise determinis- 
tic classical stochastic process in the set of all pure states 
[49, 50, 51]. The continuous time evolution is governed by an 
effective Hamiltonian, interrupted by discontinuous 'jumps' 
reflecting photon detection. This continuous part is described 
by the Schrodinger equation d t \ip) = —iH\ip) for the un- 
normalised state vector with the non-Hermitian, effective 
Hamiltonian 

H = u e \e) (e| + {lo c - ^^/2) C t c + Q(c\e) (g\ + c*\g) (e|), 

(14) 

where c is an annihilation operator for the cavity mode. The 
decreasing vector norm of \ip(t)) due to the non-unitary evo- 
lution, U (t) = e~ %m , leads to the cumulative density func- 
tion for the time at which the photodetector registers a photo- 
count, 

p(i) = i-||l#)>ll 2 - d5) 

This in turn governs the waiting time distribution in the 
stochastic process. Correspondingly, a detector click corre- 
sponds to a 'jump' in the state of the system according to 
i — > 7c|V>), where 7 = T 1 / 2 [49]. 

For the system of four atoms in cavities, with four detec- 
tors following a beam splitter network, the only change to this 
prescription is that U(t) = J|j Uj(i) and a click in detector 
k effects a 'jump' in the state according to \ty) 1 — ► dk\^), 
where dk is related to the cavity mode operators according to, 
dk = Ylj fik,jljCj> where j3 = {Pk.j) is the unitary induced 
by the beam splitter network. Note that we use j = 1 , . . . , 4 to 
label subsystems and k = 1, 4 to label detectors. 

We examine the effect of errors on the basic 4BS network 
which creates 4-GHZ states from two EPR pairs, so that the 
initial state vector is given by |^(0)) = |EPRi i2 ) |EPR 3j4 ). 
We treat the initial excitation of the protocol, |1) — > |e), as 
instantaneous and ideal, and examine the effect of mismatched 
system parameters Tj, Aj and Clj on the subsequent emission 
and detection process. The attraction of our protocol is that 
the initial state is in the two-excitation subspace, so we expect 
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6Q 2 sr 2 

FIG. 4: Contours of fidelity, F, as a function of fii = 1 + 5Qi 
(vertical axes) and fi 2 = 1 + SQ 2 (left), and r 2 = 4 + SF 2 (right). 
All other parameters are fixed at Aj — 0, ftj — 1 and Tj = 4. The 
innermost contour is at F = 1 — 1CP 4 whilst all others are at 10 -3 
intervals. 



to register exactly two detector counts. The state vector of the 
system at the end of the protocol is conditional upon which 
detectors clicked, k\ and k 2 , and at what times, t\ and t 2 — 
h + At, 

|*(ti,*i;*2,fe)) = dfc a ^(At)d fcl ^(ti)|*(0)). (16) 

This particular outcome occurs with a probability density 
function given by p(t 1 ,k 1 ;t 2 ,k2) = \\ h; h, k 2 ))\\ 2 - 
Integrating p over t and At yields the probabilities in Table I, 
as required. Note that for the ideal case, the distribution of p 
does not actually depend on fc x and fc 2 ; the only dependence 
is an overall multiplicative factor such that Table I is satisfied. 
An example is shown in Fig. 3. 

For a given combination of detectors and times, we cal- 
culate the fidelity of the resulting state with respect to 
the ideal outcomes, shown in Table I, which depend only 
on the detectors, and not the times: /(ii, ki;t2, k 2 ) = 
\{^(ki,k2)\^{ti,ki;t2,k 2 ))\ 2 (note that we have renor- 
malised the state, denoted by the lack of a tilde). In order 
to give a fair estimate of the expected fidelity, we compute the 
time-averaged fidelity for outcomes where different detectors 
click, 




where /(ti.fcijta.Aa) = *a)l*(ti, fe? *a, k 2 ))\ 2 , and 

the factor of 2 accounts for the fact that in the ideal case the 
probability of fusing the states is 1/2. / consists of a sum of 
exponentially decaying terms, so we compute F analytically, 
though the expression is rather lengthy. 

When fij = Q and Tj = V, we find F = 1, so that the 
protocol works perfectly. Otherwise, the process works with 
reduced fidelity, as shown in Fig. 4. Since we are consid- 
ering three error parameters per subsystem, we cannot show 
the dependence of F on all of them graphically, however 
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for small perturbations to the parameters, we can straightfor- 
wardly compute the dependence to quadratic order. In what 
follows, we work in units for which = 1. We find that 

F » 1 - Ei tfM.ej - Ej^ h *£M x e ja , where 



(18) 



is the vector of parametric errors in subsystem j and M s _ x are 
the coefficient matrices. For a critically damped atom-cavity 
system, T = 40, the coefficient matrices are 



128 









3 


3 


:!2 


16 


3 


9 


16 


16 . 



3 
128 












1 


1 


32 


16 


1 


3 


16 


16 



(19) 



B. Photon Loss 



Given an ideal apparatus, without any photon loss or de- 
tector failure, our fusion process would merely require four 
simple non-photon-number-resolving detectors. However, in 
practice any near-future experiment will certainly suffer sig- 
nificant photon loss. This appears potentially very damag- 
ing to our scheme (and to that of Ref. [25], but not to that 
of Ref. [24]), because we may misinterpret a photon loss as 
two photons entering a single detector. In order to counter 
this issue, we would require a limited degree of photon res- 
olution at the detectors - specifically, we must differentiate 
the three cases: 0, 1 and more than one photons. This suf- 
fices to detect a photon loss event - we would then reset the 
associated matter qubits and rebuild the graph section (analo- 
gously to Ref. [24]). Importantly, graph state fidelity will not 
be affected by undercounts, reflecting a non-unit detection ef- 
ficiency, which is the dominant problem in real world photon 
detector technologies. Undercounting is equivalent to a lossy 
channel followed by a perfect photon counter, and therefore 
detector inefficiency is simply subsumed into the total pho- 
ton loss rate. Detector over-counting, i.e. dark counts, are 
potentially harmful in the present scheme and those of Ref. 
[24, 25]. Fortunately, since we know that correct operation of 
the scheme generates precisely two photons, we will success- 
fully identify any photon loss event unless the loss occurs at 
the same time that a detector is independently subjected to a 
dark count, a process that is expected to happen with a very 
small probability. 

One could construct an adequate detector simply from two 
non-photon-resolving detectors, together with a fast switch. 
This would exploit that fact that when two photons are inci- 
dent on one detector, they are typically separated by a time 
interval of order 1/T, see Fig. 3; this time can be made long 
enough to trigger a pockels cell to redirect a possible second 
photon into a second detector, obviating the need for a true 
number-resolving detector. On occasions when the two pho- 
tons occur too close together for the second to be redirected, 
we simply undercount and assume photon loss has occurred. 

To summarize, the present scheme is potentially more sus- 
ceptible to photon loss than the 'double heralding' scheme of 



Ref. [24]. However, the issue can be dealt with using detector 
technologies that remain relatively simple - we do not require 
high fidelity photon number resolving detectors in order to 
generate high fidelity graph states. 



C. Other Errors 

This shows that the protocol is most sensitive to errors in 
the atom-cavity coupling rate, and less sensitive to detuning 
or the cavity leakage rate (the same hierarchy as observed in 
Ref. [24]) . The method used here can be adapted to include 
dephasing errors, as it was in Ref. [40], however it is rather 
more cumbersome, so for brevity we do not analyse it in detail 
here. In Ref. [40], it was found that dephasing was minimised 
when r w O, since such a critically damped system has the 
shortest lifetime of excitations in the system. It was also noted 
that dephasing was negligible when O and T are much larger 
than the dephasing rate. We expect these statements to hold 
true in this system as well, since the underlying physical pro- 
cesses are the same. 

The issue of interferometric instability is relevant to any 
scheme in which terms in the matter qubit superposition be- 
come coupled to the presence/absence of a photon in a given 
channel. Any phase noise suffered by a photon in transit 
through the apparatus ultimately can be mapped onto the mat- 
ter qubits. Fortunately, there has been enormous progress 
recently in the development of experimental techniques for 
phase locking, which should prove to be beneficial for a 
scheme of the proposed type [52]. Imperfect optical excita- 
tion can be dealt with by noting that this simply reduces the 
initial state fidelity, which thus reduces the protocol fidelity 
by an equal amount. 



VII. SUMMARY 

We have described a scheme that unifies some of the de- 
sirable features of previous work on matter qubits and graph 
states. It is able to achieve deterministic growth while using 
simple static linear optics and a 'one shot' excitation. More- 
over, the presented scheme obviates the need for continual lo- 
cal operations on qubits during graph growth, which implies 
a dramatic speedup in many systems. The scheme proves to 
have properties that make it ideal for creating the most re- 
source efficient form of graph state, the minimal graph state. 
These minimal graph states, which form the essential resource 
for a given quantum computation, without its classically effi- 
ciently tracktable Clifford-part, typically correspond to graphs 
with a high maximal vertex degree. For the preparation of 
such graph states this scheme is particularly suitable. We ob- 
serve that the use of minimal graph states is completely com- 
patible with the idea of 'just in time' entanglement generation. 
Our protocol is relatively robust to mismatch in the subsys- 
tems, and an accuracy of greater than 1% in the parameters 
will provide a fidelity of around 0.9999 in the final state. We 
hope that the scheme presented in this work can contribute to 
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bringing full-scale graph state quantum computation closer to 
practical realization. 
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